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Spatial intermittency in fully developed turbulence is an established fea- 
ture of astrophysical plasma fluctuations and in particular apparent in 
the interplanetary medium by in situ observations. In this situation the 
classical Boltzmann-Gibbs extensive thermo-statistics, applicable when 
microscopic interactions and memory are short ranged and the environ- 
ment is a continuous and differentiable manifold, fails. Upon generaliza- 
tion of the entropy function to nonextensivity, accounting for long-range 
interactions and thus for correlations in the system, it is demonstrated 
that the corresponding probability distributions (PDFs) are members 
of a family of specific power-law distributions. In particular, the result- 
ing theoretical bi-kappa functional reproduces accurately the observed 
global leptokurtic, non-Gaussian shape of the increment PDFs of charac- 
teristic solar wind variables on all scales, where nonlocality in turbulence 
is controlled via a multi-scale coupling parameter. Gradual decoupling 
is obtained by enhancing the spatial separation scale corresponding to 
increasing kappa-values in case of slow solar wind conditions where a 
Gaussian is approached in the limit of large scales. Contrary, the scaling 
properties in the high speed solar wind are predominantly governed by 
the mean energy or variance of the distribution, appearing as second pa- 
rameter in the theory. The PDFs of solar wind scalar field differences are 
computed from WIND and ACE data for different time-lags and bulk 
speeds and analyzed within the nonextensive theory, where also a partic- 
ular nonlinear dependence of the coupling parameter and variance with 
scale arises for best fitting theoretical PDFs. Consequently, nonlocality 
in fluctuations, related to both, turbulence and its large scale driving, 
should be related to long-range interactions in the context of nonexten- 
sive entropy generalization, providing fundamentally the physical back- 
ground of the observed scale dependence of fluctuations in intermittent 
space plasmas. 
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1. Introduction 



Leptokurtic, long-tailed probability distribution functions (PDF's) subject 
to a non-Gaussian core and pronounced halo are a persistent feature in a 
variety of different astrophysical environments. Those include the thermo- 
statistical properties of the interplanetary medium where the electron, 
proton and even heavy ion velocity space distributions show ubiquitousl 
suprathermal halo patterns (see Mendifffor a general review or Leubne: 
and references therein), well described by the empirical family of kappa- 
distributions, a power law in particle speed, recognized first by Vas yliun affl 
In continuation, significant progress was provided by Treumanr]5El who 
developed a kinetic theory, demonstrating that power-law velocity space 
distributions are a particular thermodynamic equilibrium state. Similarly, 
scale invariant power-law distributions are manifest in any systems rely- 
ing on self org aniz ed criticality (SOC) (BeM, Bak et alP, Watkins et alP, 
Chapman etalP3 Chapman and WatkinA^ or gravitationally bound as- 
trophysical stellar systems (Nakamichi et al!^, Chavanis and BouchelP^). 
Moreover, recently Leubner ^31 developed a theory representing accurately 
the hot plasma and dark matter density profiles in galaxies and clusters 
in the context of scale invariant power-law distributions. In all cases the 
standard Boltzmann-Gibbs-Shannon (BGS) statistics does not apply . 

Remarkably, we have to add to this diversity also the PDF's of the tur- 
bulent fluctuations of the magnetic field strength, densit y and v elocity fields 
in space and astrophysical plasmas (Leubner and Voroii^^) . In particu- 
lar, the analysis of the PDFs of the solar wind plasma is of considerable 
interest to study intermittency and multi-scale statistical properties in fully 
developed turbulence, since high resolution in situ observations are avail- 
able. The characteristics of the spectral properties of fluctuations in the 
incompressible interplanetary medium were provided in classical statistical 
theory via the phase space di stri bution obtained from ideal MHD invariants 
by Matthaeus and GoldsteirU^I and followed by a co nfirma tion of the ex- 
istence of solar wind multifractal structures (Burlag . Furthermore, 
solar wind observations were al so a ble to study the differences between fluid 
and MHD turbulence (Carbon 

The non-Gaussianity of the PDFs of the magnetic field and plasma 
fluctuations was analyzed and linked to intermittency and the fractal scal- 
ing of the solar wind MHD fluctuations (Marsch and T\|2E!3) ; followed by 
detailed investigations of the non-Gaussian fractal and/or multifractal char- 
acteristics of solar wind and related magnetospheric parameters (HnaJ2Sl). 
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Their multiscalc coupling properties and significance in view of magneto- 
spheric response was analized (Voros et alP^l Voros and JankovicovJ^j), 
including studies of multi-scale intermittency and anisotropy effects in the 
near-Earth magnetotail dynamics (Voros et alP^^j). 

WIND, ACE and Voyager observations of solar wind multi-scale statis- 
tical properties verify that the leptokurtic, long-tailed shapes of the PDFs 
at small scales represent the characteristics of intermittent turbulence and 
approach a Gaussian, reflecting a decoupled stat e, on large scales (Sorriso- 
Valvo et alPSl Burlaga et alPSl Burlaga et alPSI). In other words, the prob- 
ability of rare events is raised on small scales, where the spatial separation 
scale is characterized commonly by the differences SX(t) — X(t+r) — X(t), 
X(t) denoting any characteristic solar wind variable at time t and r is the 
time lag. Recently, intermittency was considered to appear as result of an 
interplay between stochastic Alfvenic fluctuations and coherent 2-D struc- 
tures (Bruno et al.^J) 

The empirical Castaing model (Castaing et alP3, Castaing and 
DubrullJSI, CastaingSl) introduces intermittency through fluctuations of 
log-normal distributed variances based on the idea that for constant energy 
transfer between spatial scales all variables obey a Gaussian distribution 
of fluctuations SX and hence assuming that fluctuations on different scales 
are independent. This convolution of Gaussians of different variances was 
introduced to model the non-Gaussian energy casca de character of inter- 
mittency in turbulent flows (Consolini and MicheliJ^, Guowei et alP^j 
Sorriso-ValvcE^I Schmitt and MarsarP^) where the log-normal distribu- 
tion of variances through the inertial scales was found to provide excellent 
fits to the observed leptokurtic PDFs in solar wind flows. Due to the fit- 
ting accuracy the Castaing model achieved high popularity, but appears to 
be subject to two significant shortcomings on physical grounds: the model 
provides (1) no link to non-locality and long-range interactions present in 
turbulence and (2) no justification for direct energy coupling between sep- 
arated scales. Other known shortcomings of Castaing or log-normal models 
are related to the predicted features of high-order moments , w hich violate 
the fundamental assumptions of turbulence theory (Frisch- , Arratia et 
alR. 

The global leptokurtic non-Gaussian shape of the increment PDFs re- 
quires theoretically a corresponding unique global distribution function. 
This condition can be formulated on a general level by considering the basic 
feature of turbulent flows, i.e. multi-scale coupling or nonlocality in phys- 
ical or in Fourier space, where nonlocality appears due to the presence of 
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long-range forces implying direct nonlocal interactions between large scales 
and small scales. Due to long-range interactions small and large scales are 
strongly coupled indicating that small-scale fluctuations in each time/space 
point depend on the la rge scale motions in the whole time/space domain 
and vice versa (TsinobeF^f) . Accounting for long-range interactions is a par- 
ticular feature of nonextensive systems and available from pseudo-additive 
entropy generalization. 

The classical Boltzmann-Gibbs extensive thermo-statistics constitutes a 
powerful tool when microscopic interactions and memory are short ranged 
and the environment is an Euclidean space-time, a continuous and diffcr- 
cntiable manifold. However, in the present situation we are dealing with 
astrophysical systems, generally subject to spatial or temporal long-range 
interactions evolving in a non-Euclidean, for instance multi-fractal space- 
time that makes their behavior nonextensive. A suitable generalization 
of the Boltzmann-Gibbs-Shannon entropy for statistic al e quilibrium was 
first proposed by Renyl2Hl and subsequently by TsallisSHI preserving the 
usual properties of positivity, equiprobability and irreversibility, but suit- 
ably extending the standard extensivity or additivity to nonextensivity. 
The main theorems of the classical Maxwell-Boltzmann statistics admit 
profound generalizations within nonextensive statistics (sometimes referred 
to as q-statistics where q characterizes the degree of nonextensivity of the 
system), wherefore a variety of subsequent analyses were devoted to clarify 
the mathematical and physical consequences of pseudo-additivity, for an 
early review see e.g. TsallisPH Those include a reformulation of the clas- 
sical N-body problem within the extended statistical mechanics (Plastino 
et alP) and the development of nonextensive distributions (Silva et al.™^, 
AlmeidJSI) where a deterministic connection between the generalized en- 
tropy and the resulting power-law functionals (Andrade et all^j), as well 
as the duality of nonextensive statistics were recognized (Karlin et alP^). 

Rela ting the parameters q and n by the transformation « = 1/(1 — q) 
(Leubnep^l) provided the missing link between nonextensive distributions 
and K-functions favored in space plasma physics, leading to the required 
theoretical justification for the use of ^-distributions from fundamental 
physics. Since the parameter k, a measure of the degree of nonextensivity 
of the system, is not restricted to positive values in the nonextensive con- 
text, the commonly observed core-halo twin character of the interplanetary 
electron and ion velocity space distributions was verified theoretically upon 
generalization to a bi-kappa distribution, subject to a less pronou nced c ore 
along with extended tails, as compared to a Maxwellian (LeubneiEBEHl) . 
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Recently, the PDF of the Tsallis ensemble was linked to the analysis of 
fully developed turbulence providing a relation between the nonextensive 
parameter q and the intermittency exponent m that is, a manifestation 
of multifractality of the distribution of eddies (Arimitsu and Arimits 
Arimitsu and ArimitsJ^J) as well as of scaling of the velocity structure 
functions (Arimitsu and ArimitsvJ^). Moreover, the context of general- 
ized thermo-statistics provides analytical formulas for PDFs of distance 
dependent velocity differences, linking the entropic index to the cascade 
like structure of the turbulent dynamics (BecH^). We relate in the fol- 
lowing nonlocality in turbulent flows to the presence of long-range forces 
in nonextensive systems and demonstrate in the context of entropy gener- 
alization the consistency of the the oretically derived bi-kappa distribution 
(LeubneJSHI Leubner and Voroii^J) with observed, scale dependent PDFs 
of characteristic variables in the intermittent, turbulent solar wind, where 
both, slow and high speed conditions are analyzed separately. 

2. Theory 

The standard BGS statistics is based on the extensive entropy measure 

S B = -k B },PjhiPi (1) 

where pi is the probability of the i th microstate, ks is Boltzmann's con- 
stant and Sb is extremized for equiprobability. As physical background one 
assumes that particles move independently from each other, i.e. there are 
no correlations present in the system considered. This implies isotropy of 
the velocity directions and thus the entropy appears as additive quantity 
yielding the standard Maxwellian distribution function. In other words, mi- 
croscopic interactions are short ranged and we are dealing with an Euclidean 
space time. The assumptions behind standard BGS statistics are not ap- 
plicable if one needs to account for nonlocality and long-range interactions 
in a fractal/multifractal physical environment. It is required to introduce 
correlation within the system, which is done conveniently in the context of 
nonextensive entropy generalization leading to scale-free power-law PDFs. 

Considering two sub-systems A and B one can illuminate nonextensivity 
by the property of pseudo-additivity of the entropy such that 



S K (A + B) = S K (A) + S K (B) + -S K (A)S K (B) 



(2) 
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where the entropic index k is a parameter quantifying the degree of 
nonextensivity in the system. For k = oo the last term on the right hand 
side cancels leaving the additive entropy of standard BGS statistics. Hence, 
nonlocality or long-range interactions are introduced by the multiplicative 
term accounting for correlations between the subsystems. As a measure for 
entropy mixing the entropic index n quantifies the degree of nonextensivity 
in the system and thus accounts for nonlocality and long-range interactions 
or coupling and correlations, respectively. In general, the pseudo-additive, 
K-weighted term may assume positive or negative definite values indicat- 
ing a nonextensive entropy bifurcation. Obviousely, nonextensive systems 
are subject to a dual nature since positive K-values imply the tendency to 
less organized states where the entropy increases whereas negativ n- values 
provide a higher organized state of decreased entropy, see LeubneJi^l 

The general nonextensive entropy consistent with Eq. replacing the 
classical BGS-statistics for systems subject to long-range interactions, takes 
the form (TsalliJSOl LeubneiPD) 



S K = K fc s (X>-~ 1/K -l) (3) 

In order to link the K-notation defined within — oo < k < +oo, com- 
monly applied in space plasma modeling in terms of the family of re- 
distributions, to the Tsallis q-statistics one may perform the transformation 
1/(1 - q) = k to Eq. © (Leubne]E3). k = oo corresponds to q = 1 and 
represents the extensive limit of statistical independency. Consequently, the 
interaction term in Eq. JSJl cancels recovering with respect to Eq. the 
classical Boltzmann-Gibbs-Shannon entropy Eq. Q. Eq. © applies to sys- 
tems subject to long-range interactions or memory and systems evolving 
in a non-Euclidean and multifractal space-time. A further generalization of 
Eq. |J2J for complex systems, composed of an arbitrary number of mutually 
correlated systems, is provided by Milovanov and ZelenyP^I where appro- 
priate higher order terms in the entropy appear. Once the entropy is known 
the corresponding probability distributions are available. 

In Maxwells derivation the velocity components of the distribution f(v) 
are uncorrelated where Inf can be expressed as a sum of the logarithms 
of the one dimensional distribution functions. In nonextensive systems one 
needs to introduce correlations between the components accounting for the 
long-range interactions, which is conveniently done by extremizing the en- 
tropy under conservation of mass and energy yielding the corresponding 
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one-dimensional power-law distributions as 



A 3 



1 v 2 

KVf 



(4) 



where v t corresponds to the mean energy or thermal speed of the distri- 
bution. Hence, the exponential probability function of the Maxwellian gas 
of an uncorrelated ensemble of particles is replaced by the characteristics 
of a power law where the sign of k, indicated by superscripts, governs the 
corresponding entropy bifurcation. Note that the distribution @ can be 
derived entirely using general methods of statistical mechanics without in- 
troducing a specific form for long-range interactions. Incorporating the sign 
of k into Eq. generates a dual solution with regard to positive and neg- 
ative K-values, respectively, resulting also in two different normalizations 
A ± . 

The entropy bifurcation appears to be manifested also in higher order 
moments yielding f or th e second moments K-dependent generalized temper- 
atures, see LeubneiEHl Furthermore, the positive solution is restricted to 
k > 3/2 whereas the negative solutions are subject to a cutoff in the distri- 
bution at v max = vt\f[n). Both functions, / + and /~ in Eqs. Q approach 
the same Maxwellian as k — ► oo. Fig. 1, left panel, demonstrates schemati- 
cally the non-thermal behavior of both, the suprathermal halo component 
and the reduced core distribution, subject to finite support in velocity space, 
where the case k = oo recovers the Maxwellian equilibrium distribution. 

Any unique and physically relevant nonextensive PDF must obey the 
following three conditions: (a) the distribution approaches one and the same 
Maxwellian as k — > oo, (b) a unique, global distribution must be definable 
by one single density and a unique temperature and (c) upon variation of 
the coupling parameter k particle conservation and adiabatic evolution are 
required, such that a redistribution in a box (a source free environment) 
can be performed. Subject to these constraints the appropriate mathemat- 
ical functional, representing observed core-halo structures in nonextensive 
astrophysical environments, is available from the elementary combination 
f c h = B c h(fh + fc), B c h being a proper normalization constant. In this con- 
text the full velocity space bi-kappa distribution, compatible with nonex- 
tensive entropy generalization and obeying the above constraints, reads 



F ch (v;n) 



N 



v t 



-G{k) 



lv 2 



lv 2 



(5) 
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2 2 
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Fig. 1. Left panel: A schematic plot of the characteristics of the nonextensive bi-kappa 
distribution family: with k = 3 the outermost and innermost curve correspond to the halo 
fh and core f c distribution fraction. For increasing K-values both sets of curves merge at 
the same Maxwellian limit, indicated as bold line, from outside and f c from inside. 
Right panel: A nonextensive bi-kappa fit (solid line) with n = 1.8 of an observed PDF 
(dashed line) obtained from ACE magnetic field amplitude data. A Gaussian (dashed- 
dotted line) and a conventional n— distribution are provided for comparison. 



The last term on the right-hand side denotes an expression subject to 
a thermal cutoff at the maximum allowed velocity v max = \fR,v t , which 
limits also required integrations. Fig. 1 (right panel) indicates also that the 
core of the distribution, corresponding to the last term on the right-hand 
side of Eq. (JSJ, contributes only for small SB. For details regarding the 
correspon ding second moments or generalized temperatures see LeubneJ3Hl 
or LeubnerS21 The function G(k) is defined by 



G{k) = 



K 1 ' 2 T(k- 1/2) k 1 ' 2 T{k+1) 



I» 



r(/s + 3/2) 



(6) 



from the normalization and is subject to a particular weak K-dependence 
where G(k) ~ 1/2, see Leubnei@Slfor a graphical illustration and discussion. 
Hence, the normalization is independent of the parameter k and the factor 
1/2 reflects consistently the superposition of the two counter-organizing 
contributions in Eq. J^J). For k — oo, G(k) — 1/2 and the power laws 
in the brackets of the right hand side of Eq. © turn each into the same 
Maxwellian exponential. 

With regard to this specific mathematical feature the energy levels Ei 
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9 



of the turbulent spectrum can be related to the corresponding kinetic en- 
ergy of velocity differences Sv(t) = v(t + r) — v(t) between two points of 
separation r allowing to transform the ID Maxwellian particle distributin 
of mean energy v t into the mathematical form of a Gaussian of variance a. 
Upon normalizing the one-dimensional bi-kappa particle distribution (JSJ) to 
unity and assigning the distribution variance a to the thermal spread Vt the 
'Maxwellian form' of the bi-kappa distribution transforms to a 'Gaussian 
form' of a global bi-kappa PDF as 



This two-parameter PDF is applicable to the differences of the fluctua- 
tions 5X(t) = X(t+r) — X(t) of any physical variable X in the astrophysical 
system considered. Here k assumes a physical interpretation defining the 
degree of nonextensivity or nonlocality in the system, thus being a mea- 
sure of th e d egree of organization or intermittency, respectively (Leubner 
and VoroJISI) and a denotes the distribution variance. Again, depending 
on the particular choice of the parameters, the last term on the right-hand 
side is subject to a cutoff, treated by limiting the integration to the proper 
interval. Fig. 1, left panel, illuminates also that large values of 5X, corre- 
sponding to the tails of the distribution, are represented by the first term 
on the right-hand side of Equation J7J). As n — > oo the bi-kappa distribution 
P c h{SX; n,a) approaches a single Gaussian. 

Note, that Equation 5 represents the full velocity space bi-kappa distri- 
bution, while Equation 7 is related to the time shifted differences of any 
physical variale X. Physically, Equation 5 describes particles which motion 
is controlled by long-range forces and interactions. The basic assumption for 
deriving the velocity space bi-kappa distribution was the pseudo-additivity 
of the entropy of particle sub-systems expressed through Equations 2 and 
3. It is important to recognize that, the same type of expression for bi- 
kappa distribution is obtained, if instead of interacting particles we assume 
interacting coherent structures with the same pseudo-additivity property 
of the entropy as in Equation 2. In the context of MHD, nonpropagating 
multi-scale coherent structures or flux tubes can interact, deform and pro- 
duce new sites of nonpropagating fluctuations. Coherent structures of the 
same polarity merge into a structure with lower local energetic stat e, while 
structures of opposite polarities may repel each other (ChangSS Bruno 
and Carbone^). These coherent structures can be considered as discrete 
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interacting "particles" in MHD flows, responsible for the particular entropy 
within the system and validating the analogy to the kinetic level of PDFs 
(Leubner and Voros^j) . A usual way of statistical analysis of intermittency 
due to the occurrence of coherent structures in turbulence is through two- 
point differences of fluctuations. Therefore, in the expression of bi-kappa 
distribution for turbulent interactions (Equation 7) differences of physical 
variables can be used instead of the velocity notation. Moreover, passive 
scalars as the magnetic field follow the dynamics of V or SV (see Voros et 
alP). 

One of the basic future of turbulent flows is multi-scale redistribution 
of energy during which interacting coherent structures appear, reducing 
the entropy of the system and leading also to negative kappa values. At 
the same time, turbulence enhances dissipation and mixing of the plasma, 
which increases entropy and can be described in terms of positive kappa 
values. The presence of processes which increase entropy and those which 
decrease entropy in turbulence advocates therefore the introduction of bi- 
kappa distributions. 

We proof in the following the relevance of the nonextensive, global bi- 
kappa PDF J7J) on the observed scale dependence of the PDFs of the differ- 
ences of magnetic field, velocity and density variables in the intermittent, 
turbulent interplanetary medium. 

3. Application: Scale dependent interplanetary PDFs 

Based on the nonextensive two parameter bi-kappa distribution J7J we com- 
pare the PDFs obtained from slow and fast solar wind data with particular 
attention to the scale dependent changes of the two physically interpretable 
parameters (k, a) involved. We do not consider any occurrence of discon- 
tinuities or shocks in the system. The problem of interaction of turbulence 
with large scale structures as shocks is investigated elsewhere. (Voros et 
alP) 

For each data set the magnetic field and plasma parameter increments 
were calculated at a given time lag r by SX(t) — X(t + r) — X(t), where 
X represents any physical variable considered. For each realization the em- 
pirical probability distribution function (histogram) was then computed. 
SX(t) is binned into n equal spaced boxes and the number of elements in 
each box was computed where the robustness of the histograms against n 
is tested. SX(t) represents characteristic fluctuations at the time scale r 
or, equivalently, across eddies of size I ~ vt. Hence, by changing r it is 
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Fig. 2. The PDF of the increments of observed ACE magnetic field fluctuations for 
r = 100 and a resolution of 16s as compared to the bi-kappa function with re = 1.8. Based 
on the same data the central panel provides the characteristics for increased t = 2000 
where k assumes a value of 3.0 for the best representation. The PDF of large-scale 
magnetic field fluctuations, r = 10000, are well modeled by a Gaussian with n = oo, 
right panel. The dotted lines correspond to the standard deviations. 

possible to analyze the statistical features of fluctuations in different time 
scales, which roughly correspond to those statistical characteristics across 
turbulent eddies of size I ~ vt. 

In the following we demonstrate by means of Eq. (J7J) from first principle 
statistics that the strong non-Gaussianity of the PDF of small scale fluc- 
tuations should be associated physically with long-range interactions pro- 
vided in nonextensive systems by pseudo-additive entropy generalization. 
The scale dependence of the PDF in the solar wind can be represented 
accurately via the tuning parameters k and a of the bi-kappa functional. 

In Fig. 1, right panel, we focus on the bi-kappa model demonstrat- 
ing that the nonextensive context, generates a precise representation for 
the observed solar wind PDF characterizing the intermittency of the small 
scale fluctuations. For comparison also a Gaussian and the conventional 
k — function (LeubneJ^) , subject to the same k— value but not able to re- 
produce the structure of the PDF of small scale fluctuations, are provided. 
It is also possible to assume that the smallest fluctuations are random and 
uncorrelated, and this might be the reason why the central part of the distri- 
bution near the maximum is well-fitted by a Gaussian, as in the right panel 
in Fig. 1. Note, however, that the characteristic scale of the differences of 
fluctuations is introduced through the time lag r. Therefore, SB — > near 
the maximum of the distribution simply means that, at the scale r the 
differences between the corresponding values of B(t + r) and B(t) are very 
small. Combined with the finite precision of PDFs estimation, small two- 
point differences can really produce uncorrelated fluctuations near PDFs 
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Fig. 3. The PDF of the increments of observed WIND velocity field magnitude fluc- 
tuations for t = 10 and a resolution of 92s as compared to the bi-kappa function with 
k = 2. Based on the same data the central panel provides the characteristics for in- 
creased t — 70 where n assumes a value of 3.5 for the best representation. The PDF 
of large-scale velocity fluctuations are well modeled by a Gaussian with k = oo, right 
panel. The dotted lines correspond to the standard deviations. 

maxima. This does not mean, however, that the Gaussian distribution is 
the correct answer for the proper description of the central part of the 
distribution at the small scales r. If it was true, we should suppose that, 
at smaller and smaller scales r, the fluctuations become more and more 
random and we are getting closer and closer to a Gaussian distribution. 
Actually the opposite is true, the peakedness of the distribution increases 
as the scale r decreases (see later). 

3.1. Slow speed solar wind 

As an example, in Fig. 2 undisturbed solar wind ACE magnetic field am- 
plitude data of 16s time resolution are analyzed where the dimcnsionlcss 
t is multiplied by the resolution to generate an effective time-lag. In par- 
ticular, the scale dependent PDF evolution of magnetic field fluctuations 
is subject to a two point separation scale of t = 100, 2000 and 10000. 
The corresponding best fits of the bi-kappa distribution are obtained for 
k = 1.8,3 and oo, measuring the degree of nonextensivity, or coupling, re- 
spectively, through long-range interactions and the dotted lines refer to the 
standard deviation. The accuracy of the bi-kappa distribution fit demon- 
strates that non-locality in turbulence, when introduced theoretically by 
long-range interactions through the nonextensive context, generates a pre- 
cise representation for the observed PDFs characterizing the intermittency 
of the fluctuations at all scales. 

Based on WIND velocity field magnitude data Fig. 3 presents an anal- 
ysis of the scale dependence of interplanetary PDFs of the velocity field 
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Fig. 4. Left panel: The PDF of the increments of observed WIND density fluctuations 
for t = 10 and a resolution of 92 sec. as compared to the bi-kappa function with k = 2. 
Based on the same data the central panel provides the characteristics for increase r = 70 
where k assumes a value of 3.5 for the best representation. The PDF of large scale density 
fluctuations with r = 900 are well modeled by a Gaussian with k = oo, right panel. 



magnitude showing in three plots from left to right the decreasing kurto- 
sis with increasing time-lags, where the observational uncertainty is again 
indicated by the standard deviation (dotted lines). Constraint by a time 
resolution of 92sec the two point time separation r assumes the values 
10, 70 and 900 from left to right. The solid lines represent best fits to the 
observed PDFs with nonextensive distributions, where the corresponding 
k is determined as k — 2,3.5 and oo. The strong non-Gaussian character 
of the leptokurtic PDFs (left panel) , exhibiting pronounced tails associated 
with solar wind turbulence and intermittency in small-scale fluctuations, 
finds again an accurate analytical fit and hence a physical background in 
the nonextensive representation. The non-Gaussian structure is somewhat 
softened for enhanced t — 70 (central panel) but again precisely mod- 
eled within the pseudo-additive entropy context, turning into the Gaussian 
shape of large scale fluctuations, which is independent of the increment 
field. 

Fig. 4 provides the corresponding nonextensive analysis of the scale 
dependence of the density fluctuations obtained from WIND data (92 sec 
resolution) . Evidently, the scale dependent characteristics of the observed 
PDFs of the increment fields SX(t) = X(t + t) — X(t) for all solar wind 
variables evolve simultaneously on small scales approaching independency 
of the increment field in the large scale Gaussian. With separation scales of 
r = 10, 70, 900 the corresponding evolution of the PDFs of observed density 
fluctuations are best represented by the same values of n = 2, 3.5, oo, as for 
the velocity field magnitude, accounting for nonlocal interactions in the 
nonextensive theoretical approach. 
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Fig. 5. The PDF of the increments of observed ACE high-speed associated magnetic 
field magnitude fluctuations (16 s time resolution). Top-left: fluctuations at the scale 
r = 10 as compared to the bi-kappa function with k = 1.4 and a = 0.05; Top-right: 
r = 40, k = 1.4 and a = 0.12; Left-bottom: r = 400, k = 1.2 and a = 0.6; Left-right: 
t = 10000, K = oo and cr = 15 (Gaussian fit). 



3.2. High speed solar wind 

The four panels in Fig. 5 show PDFs of high speed associated magnetic field 
magnitude fluctuations. The two-point statistics is demonstrated in the 
subplots from top-left to bottom-right for the scales r = 10,40,400, 10000. 
The effective time-lag is obtained after multiplying r with time resolution of 
lQsec. The corresponding best fits reveal differing statistical features of high 
speed associated magnetic fluctuations. In comparison with low speed data 



5, 2008 11:28 



WSPC/Trim Size: 9in x 6in for Review Volume 



Leubner 



Nonextensive entropy and space plasma turbulence 



15 




-5 5 

8 N [normalized] 



■5 5 

8 N [normalized] 




-5 5 

8 N [normalized] 



■5 5 

8 N [normalized] 



Fig. 6. The PDF of the increments of observed ACE high-speed associated density 
fluctuations (64 s time resolution). Top-left: fluctuations at the scale r = 10 as compared 
to the bi-kappa function with k = 1.5 and a = 0.3; Top-right: r = 40, n = 1.5 and 
<t = 0.5; Left-bottom: r = 400, n = 1.05 and cr = 0.9; Left-right: r = 3000, k = oo and 
<r = 25 (Gaussian fit). 



the degree of nonextensivity does not change during high speed intervals, 
k = 1.4,1.4,1.2, over the range of scales r = 10,40,400, and only for 
t = 10000, k reaches oo. In contrast to the slow wind data, however, good 
quality high speed fits can be achieved only when the standard deviation a 
is changed. 

It indicates that the abundance of large scale energy content of high 
speed flows may facilitate to maintain the degree of nonextensivity and 
self-organisation unchanged over the considered scales. On the other hand, 
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we have no clear explanation yet for the observed changes of the standard 
deviation. We can speculate that changes in sigma appear because of the 
ample changes in the amplitudes of two-point fluctuations in the solar wind, 
having solar origin or being generated by local processes absent in the slow 
wind. Obviously, further comparative case studies on multi-scale fluctua- 
tions are needed to answer the question of the relative contribution of local 
processes versus processes originating in the solar corona to the observed 
behaviour of statistical moments in the fast and slow solar wind. 

The four panels in Figs. 6 and 7 provide the same qualitative behavior for 
high speed associated ACE density (64sec resolution) and WIND magnetic 
field (3sec resolution). 

4. Discussion and conclusions 

The Wind and ACE solar wind data analysis unambiguously manifests that 
the PDFs of large scale density, velocity and magnetic field fluctuations are 
well represented by a Gaussian, turning into leptokurtic peaked distribu- 
tions of strong non-Gaussianity in the center along with a pronounced tail 
structure at smaller scales. In particular, the PDFs of large-scale magnetic 
field fluctuations, not related to the increment field are known to be subject 
to relatively small deviations from the Gaussian statistics and are well fit- 
ted by the Castaing distribution, a convolution of Gaussians with variances 
distributed according to a log-normal distribution (CastainjpS Padhye et 
al!""|). Assuming a constant energy transfer rate between spatial scales all 
quantities exhibit a Gaussian distribution of fluctuations in this context. 
Independent of the physical situation considered, the Castaing distribution 
provides a multi-parameter description of observed PDFs, plausible in this 
case, since the large-scale fluctuations of the interplanetary magnetic field 
are generated by a variety of discrete coronal sources. If individual coro- 
nal sources evoke Gaussian distributed magnetic fields, the net magnetic 
fluctuations can be modeled by their superposition with a spread of the 
corresponding variances. 

Contrary, small-scale fluctuations are associated with local intermittent 
flows where fluctuations are concentrated in limited space volumes. Conse- 
quently, the PDFs are scale dependent and intermittency generates long- 
tailed distributions. It is customary to use n — th order absolute powers of 
the plasma variables and mag netic field increments (n — th order structure 
functions (Marsch and Tu 22 , Pagel and Balog allowing to investigate 
the multi-scale scaling features of fluctuations. Direct studies of observed 
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Fig. 7. The PDF of the increments of observed ACE high-speed associated magnetic 
field magnitude fluctuations (3 s time resolution). Top-left: fluctuations at the scale 
t = 10 as compared to the bi-kappa function with k = 1.4 and a = 0.003; Top-right: 
T = 40, k = 1.4 and a = 0.006; Left-bottom: r = 400, k = 1.8 and a = 0.06; Left-right: 
t = 4000, K — oo and o~ = 1 (Gaussian fit). 

PDFs of the increment fields 6X(t) = X(t + r) — X(t) for any character- 
istic solar wind variable at time t and time lag r revealed departures from 
a Gaussian distribution over multiple scales (Sorriso-Valvo»"|) and a n in - 
crease of kurtosis (intermittency) towards small scales (Marsch and 
The PDFs are also found to be leptokurtic, which indicates the turbulent 
character of the underlying fluctuations. Sorriso-Valvo et alP^lhave shown 
that the non-Gaussian behavior of small-scale velocity and magnetic field 
fluctuations in the solar wind can also be described well by a Castaing dis- 
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Fig. 8. The functional dependence of the spatial scale t(k) (left panel) for slow speed 
and t(ct) (right panel) for high speed conditions. 

tribution where the individual sources of Gaussian fluctuations appear at 
small-scales in turbulent cascades. 

From the corresponding nonextensive WIND data an alys is of the den- 
sity and magnetic field fluctuations (Leubner and Voroi^) it is evident 
that the scale dependent characteristics of the observed PDFs of the in- 
crement fields SX(t) = X(t + t) — X(t) for solar wind variables evolve 
simultaneously on small scales, approaching independency of the increment 
field in the large- scale Gaussian. Highly accurately, the overall scale de- 
pendence appears as a general characteristic of quiet astrophysical plasma 
environments indicating a universal scaling dependence between density, 
velocity and magnetic field intermittency within the experimental uncer- 
tainties. This strong correlation implies that the scale dependencies of all 
physical variables are coupled, where the solar wind Alfvenic fluctuations 
provide a physical basis of the velocity and magnetic field correlations. On 
the other hand, according to recent analyses, the magnetic field intensity 
exhibits a higher degree of intermittency than the s olar wind bulk veloc- 
ity, both in fast and slow winds (Sorriso-Valvo et alP^J). However, Veltri 
and Mangene}^^ found that the most intermittent structures in the slow 
wind are shock waves, displaying similar intermittency in the magnetic field 
intensity and bulk velocity. Furthermore, the proportionality between den- 
sity fluctuations and the magnetic field and velocity fluctuations is already 
maintained in the solar wind by the presence of weak spatial gradients 
(SpangleJSH). 

Fig. 8 (left panel) provides an estimation of the functional dependence 
between the time lag r and the nonextensive parameter k for best fitting 



February 5, 2008 11:28 



WSPC/Trim Size: 9in x 6in for Review Volume 



Leubner 



Nonextensive entropy and space plasma turbulence 19 

bi-kappa functions to the observed PDFs for slow speed solar wind condi- 
tions (v < 400fcms~ 1 ). As significant global behavior the scale dependence 
of the PDFs for quiet conditions appears to be independent of the variance 
or mean energy of the distribution. The best fitting bi-kappa functions are 
found for a constant corresponding parameter a, indicating that the pa- 
rameter k, measuring the degree of coupling within the system, governes 
primarily the scale dependence of the PDF in the slow solar wind. As r in- 
creases from small scales to the intermediate regime a pronounced plateau 
formation is established, i.e. the relative increase in K-values with enhanced 
scales appears reduced. In other words, the PDF shape appears at inter- 
mediate scales to be independent of the spatial separation scale. Such a 
behavior may indicate the presence of a transitional dynamical element 
characterizing a balance between long and short-range interactions. Con- 
trary, in high speed streams (v > 400fcms~ 1 ) best fits of bi-kappa functions 
to the observed scale dependent PDFs are found when keeping k constant 
and varying only the parameter a, left panel in Fig. 8. Hence the scaling 
features in the fast wind appear independent of the degree of coupling con- 
trolled by k, but rely predominantly on changes of the variance or mean 
energy. In summary, the theoretical nonextensive context indicates a sig- 
nificant and physically contrary scaling behavior in slow and fast wind. For 
slow solar wind conditions the correlations/intermittency governed by k de- 
crease with increasing scale, whereas the characteristic energy governed by 
a remains constant. Contrary the scaling properties in high speed streams 
are characterized by constant correlations (n) but enhanced variance with 
increasing scale. 

k — distributions reproduce the Maxwell-Boltzmann distribution for k — > 
oo, a situation identifying k as an ordering parameter that acounts for cor- 
relations within the system. Highly correlated turbulent conditions charac- 
terized by kappa distributions represent stationary states far from equilib- 
rium where a generalization of the Boltzmann-Shannon entropy, as mea- 
sure of the level of organization or intermittency, applies (Goldstein and 
LebowitsPI TreumanrP^J) . Physically this can be understood considering 
a system at a certain nonlinear stage where turbulence may reach a state 
of high energy level that is balanced by turbulent dissipation. In this en- 
vironment equilibrium statistics can be extended to dissipative systems, 
approaching a stationary state beyond thermal equilibrium (Gotoh and 
Kraichnan ). Since turbulence is driven in the solar wind by velocity 
shears we have choosen for the data analysis intervals of low speed so- 
lar wind with limits in velocity space, where the driving and dissipation 
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conditions do not change significantly, maintaining therefore the dynami- 
cal equilibrium condition. High speed intervals may depend more on the 
occurence of dynamical processes on the Sun. 

A multi-scale cascade mechanism is not the only way for a realization of 
long range interactions. Let us provide a physical situation where large and 
small scales are directly coupled and the nonlocal energy transfer is no t in - 
duced by cascading processes, following e.g. a log-normal model (Friscrl^|). 
An example is found in the context of MHD where Chang^l and Con- 
solini and Chang^ proposed an intermittent turbulence model for the 
solar wind and for the Earth's magnetotail, which comprises neither cas- 
cades nor requires local interactions in Fourier space. In this scenario non- 
propagating or convected fluctuations generate multiscale coherent struc- 
tures (e.g. flux tubes), which can interact, deform and produce new sites 
of non-propagating fluctuations. Coherent structures of the same polarity 
merge into a structure with lower local energetic state, while structures of 
opposite polarities may repel each other. This coherent structures can be 
considered as discrete interacting 'particles' in MHD flows, responsible for 
the particular entropy within the syst em and validating the analogy to the 
kinetic level of PDFs. Chang et alP^have computed PDFs of the intermit- 
tent fluctuations from direct numerical simulations of interacting coherent 
structures. The resulting PDFs have typical leptokurtic shapes, which can 
be well fitted again by a variety of models, including those with predomi- 
nantly local interactions in Fourier space. Since all models provide similar 
fitting accuaracy it is required to focus on the underlying physical situation 
in turbulent flows. Hence, with regard to nonlocal interactions not based on 
cascade processes the nonextensive entropy approach provides physically a 
justification for nonlocal interactions and should therefore be favored over 
cascade models in such processes. 

The entropy quantifies the degree of structuring in intermittent tur- 
bulence expressed through singular multifractal measures, where also the 
parameter n (or q) i s related to the extrem es o f multifractal distributions 
(Lyra and TsalliJ^, Arimitsu and ArimitsrP^I Arimitsu and ArimitsvJ^, 
BeclPl). Since the nonextensive entropy approach is independent of the 
mechanism leading to the structures - in both situations, cascading pro- 
cesses and multiscale interacting coherent structures - or even in coexisting 
situations, the entropy concept can be applied for the analysis and quantifi- 
cation of resulting characteristics in turbulence. In summary, the majority 
of hitherto existing models of intermittency in the solar wind essentially 
correspond to the cascade picture of turbulence. Small-scale intermittency, 
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however, can be associated also by emerging topological complexity of co- 
herent sturctures in turbulence, which might be understood better through 
entropy concepts, disregarding the goodness criteria of different fits. 

We provided specific examples based on multiscale interacting coherent 
structures where the traditional cascade - and thus e.g. the log-normal ap- 
proach beside others - should not be applied for physical reasons (not in 
terms of PDF fitting accuracy). Therefore the proposed context based on 
entropy ganeralization has a potential to describe the underlying physics 
suitably and thus justifies the nonextensive approach. Certainly this must 
be viewed as an incomple te c oncept in view of the complexity of intermit- 
tence in turbulence (Coher|SHl) ) i n particular regarding the interactions with 
large scale structures. 

Summarizing, a bi-kappa distribution family turns out theoretically as 
consequence of the entropy generalization in nonextensive thermo-statistics. 
The two-parameter global bi-kappa function provides theoretically access 
to the scale dependence of the PDFs observed in astrophysical plasma tur- 
bulence. The redistribution of a Gaussian on large scales into highly non- 
Gaussian leptokurtic and long-tailed structures, manifest on small scales, 
is theoratically well described by the family of nonextensive distributions. 
Pseudo-additive entropy generalization provides the required physical in- 
terpretation of the parameter k in terms of the degree of nonextensivity of 
the system as a measure of nonlocality or coupling due to long-range inter- 
actions whereas the variance a measures the mean energy in the system. 
The scale dependence in the slow speed solar wind is sensitive to variations 
of k and in high speed streams to variations of a. We argue that multi- 
scale coupling and intcrmittcncy of the turbulent solar wind fluctuations 
must be related to the nonextensive character of the interplanetary medium 
accounting for long-range interaction via the entropy generalization. 
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